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Anomalous Diffusion on the Hanoi Networks 
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Abstract. - Diffusion is modeled on the recently proposed Hanoi networks by studying the mean- 
square displacement of random walks with time, (r 2 ) ~ £ 2 / d «> Jt is found that diffusion - 
the quintessential mode of transport throughout Nature - proceeds faster than ordinary, in one 
case with an exact, anomalous exponent d w = 2 — log 2 (</>) = 1.30576.... It is an instance of 
a physical exponent containing the "golden ratio" (j> = (l + v§) /2 that is intimately related 
to Fibonacci sequences and since Euclid's time has been found to be fundamental throughout 
geometry, architecture, art, and Nature itself. It originates from a singular renormalization group 
fixed point with a subtle boundary layer, for whose resolution <j> is the main protagonist. The 
origin of this rare singularity is easily understood in terms of the physics of the process. Yet, 
the connection between network geometry and the emergence of <j> in this context remains elusive. 
These results provide an accurate test of recently proposed universal scaling forms for first passage 
times. 



Introduction. — The study of anomalous diffusion 
is an integral part in the analysis of transport processes 
in complex materials [1-5]. Random environments often 
slow transport significantly, leading to sub-diffusive be- 
havior. Much attention has thus been paid to model sub- 
diffusion on designed structures with some of the trap- 
pings of disordered materials, exemplified by Refs. [6-11]. 
Even self-organized critical processes can be shown to 
evolve sub-diffusively, controlled by the memory of all 
past events [12]. On the other hand, tracer particles in 
rapidly driven fluids may exhibit super-diffusive behav- 
ior [13], typically modeled in terms of Levy flights [3, 14]. 
Both regimes are self-similar, fractal generalizations of or- 
dinary diffusion. 

In this Letter we consider diffusion on two new net- 
works, which yield interesting realizations of super- 
diffusive behavior. Both of these networks were introduced 
to explore certain aspects of small- world behavior [15]. 
Their key distinguishing characteristic is their ability to 
mix a geometric backbone, i. e. a one-dimensional lat- 
tice, with small-world links in a non-random, hierarchical 
structure. In particular, these networks permit a smooth 
interpolation between finite-dimensional and mean-field 
properties, which is absent from the renormalization group 
(RG) due to Migdal and Kadanoff, for instance [16]. The 
unusual structure of these networks recasts the RG into 



a novel form, where the equations are essentially those 
of a one-dimensional model in which the complex hierar- 
chy enters at each RG-step as a (previously unrenormal- 
ized) source term. This effect is most apparent in the 
real-space RG for the Ising models discussed in Ref. [15]. 
It is obscured in our dynamic RG treatment below, since 
these walks are always embedded on the lattice backbone. 
On the practical side, their regular, hierarchical structure 
allows for easily engineered implementations, say, to effi- 
ciently synchronize communication networks [15]. Regard- 
ing diffusion, one of the networks proves to be merely an 
incarnation of a Weierstrass random walk found for Levy 
flights [14] with ballistic transport, while the other net- 
work shows highly non-trivial transport properties, very 
much unlike a Levy flight, as revealed by our exact RG 
treatment. The fixed point equations are singular and ex- 
hibit a boundary layer [17]. It provides a tangible case of 
a singularity in the RG [10, 18] that is easily interpreted 
in terms of the physics. 

Generating Hanoi Networks. — In the Tower-of- 
Hanoi problem [19], disks of increasing size, labeled i = 1 
to k from top to bottom, are stacked up and have to be 
moved in a Sisyphean task into a 2nd stack, disk-by-disk, 
while at no time a larger disk can be placed onto a smaller 
one. To this end, a 3rd stack is provided as overflow. First, 
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Fig. 1: Depiction of the planar "Tower-of-Hanoi" network HN3. 
Here, the Id— backbone of sites extends over < n < oo. 

disk 1 moves to the overflow and disk 2 onto the 2nd stack, 
followed by disk 1 on top of 2. Now, disk 3 can move to 
the overflow, disk 1 back onto disk 4, disk 2 onto 3, and 
1 onto 2. Now we have a new stack of disks 1, 2, and 
3 in prefect order, and only k — 3 more disk to go! But 
note the values of disk-label i in the sequence of moves: 
1-2-1-3-1-2-1-4-1-2-1-3-1-2-1-5-..., and so on. 

Inspired by models of ultra-slow diffusion [6,7], we cre- 
ate our networks as follows. First, we lay out this sequence 
on a Id— line of nearest-neighbor connected sites labeled 
from n = lton = L = 2 k — 1 (the number of moves re- 
quired to finish the problem). In general, any site n(^= 0) 
can be described uniquely by 

n = 2 i - 1 (2i + l), (1) 

where i is the label of the disk moved at step n in the 
sequence above and j = 0, 1, 2 . . .. To wit, let us further 
connect each site n to the closest site n' that is 2* steps 
away and possesses the same value of i, both only having a 
site of value at most i + 1 between them. According to the 
sequence, site n = 1 (with i = 1) is now also connected to 
n' = 3, 5 to 7, 9 to 11, etc. For sites with i = 2, site n = 2 
now also connects to n' = 6, 10 to 14, 18 to 22, etc, and so 
on also for i > 2. As a result, we get the network depicted 
in Fig. 1 that we call HN3. Except at the boundary, each 
site now has three neighbors, left and right along the Id 
"backbone" and a 3rd link to a site 2* steps away. If we 
further connect each site also to a fourth site 2 l steps in 
the other direction and allow j = 0, ±1, ±2, . . ., we obtain 
the network in Fig. 2, called HN4, where each site now has 
four neighbors. 

These new "Tower-of-Hanoi" networks - a mix of local, 
geometric connections and "small- world" -like long-range 
jumps - has fascinating properties. It is recursively de- 
fined with obvious fractal features. A collection of the 
structural and dynamic features of HN3 and HN4 are dis- 
cussed in Ref. [15], such as results for Ising models and 
synchronization. 




Fig. 2: Depiction of the "Tower-of-Hanoi" network HN4. Here, 
the Id— backbone of sites extends over — oo < n < oo. The site 
n = 0, not covered by Eq. (1), is special and is connected to 
itself here. 

Diffusion on the Hanoi Networks. To model dif- 
fusion on these networks, we study simple random walks 
with nearest-neighbor jumps along the available links, but 
using the one-dimensional lattice backbone as our metric 
to measure distances, which implies a fractal dimension 
of df = 1. Embedded in that space, we want to calcu- 
late the non-trivial diffusion exponent d w defined by the 
asymptotic mean-square displacement 

(r 2 )~t 2/4 . (2) 

A more extensive treatment yielding also first-return prob- 
abilities is given elsewhere [20] . 

First, we consider a random walk on HN4. The "master- 
equation" [21] for the probability of the walker to be at 
site ro, as defined in Eq. (1), at time t is given by 

Vn,t = ^^-[Pn-l,t-l + Vn+l,t-l] 

+f[P»-2«,t-l+7W,t-l], (3) 

where p is the probability to make a long-range jump. 
(Throughout this Letter, we considered p uniform, inde- 
pendent of n or t). A detailed treatment of this equa- 
tion in terms of generating functions is quite involved and 
proved fruitless, as will be shown elsewhere [20]. Instead, 
we note that the long-time behavior is dominated by the 
long-range jumps, as discussed below for HN3. To simplify 
matters, we set p = 1/2 here, although any other finite 
probability should lead to the same result. We make an 
"annealed" approximation, i. e., we assume that we hap- 
pen to be at some site n in Eq. (1) with probability 1/2 1 , 
corresponding to the relative frequency of such a site, yet 
independent of update-time or history. This ignores the 
fact that in the network geometry a long jump of length 
2 l can be followed only by another jump of that length or 
a jump of unit length, and that many intervening steps 
are necessary to make a jump of length 2 4+1 , for instance. 
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Here, at each instant the walker jumps a distance 2 l left 
or right irrespectively with probability 1/2 4 , and we can 
write 



*Pn,t — ^ ^ r T n n >'Pn',t— 1 



(4) 



with 



T > 

J- n,n' 



a- 1 
2a 



n—n'b i +6 

n—n, — b 

0. (5) 



where a = 6 = 2. Eqs. (4-5) are identical to the Weier- 
strass random walk discussed in Refs. [14,22] for arbitrary 
1 < a < b 2 . There, it was shown that d w = ln(a)/ln(6), 
which leads to the conclusion that d w = 1 in Eq. (2) for 
HN4, as has been predicted (with logarithmic corrections) 
on the basis of numerical simulations in Rcf. [15]. These 
logarithmic corrections are typical for walks with marginal 
recurrence, which typically occurs when d w — df, such as 
for ordinary diffusion in two dimensions [4] . 

For HN3, the master-equation in the bulk reads for 

1 -p 

Pn,t = —^—[Pn-l,t-l+'P n+ l,t-l]+p'P n > y t-l, 



n' = < 



j even, 



n — 2 l , j odd, 



(6) 



with n as in Eq. (1), and p as before. 

In the RG [21,23] solution of Eq. (6), at each step we 
eliminate all odd sites, i. e., those sites with i = in 
Eq. (1). As shown in Fig. 3, the elementary unit of sites 
effected is centered at all sites n having i = 1 in Eq. (1). 
We know that such a site n is surrounded by two sites of 
odd index, which are mutually linked. Furthermore, n is 
linked by a long-distance jump to a site also of type i = 1 
at n ± 4 in the neighboring elementary unit, where the 
direction does not matter here. The sites n±2, which are 
shared at the boundary between such neighboring units 
also have even index, but their value of i > 2 is indeter- 
mined and irrelevant for the immediate RG step, as they 
have a long-distance jump to some sites m± at least eight 
sites away. 

Using a standard generating function [21], 



x„(z) 



OG 

t=0 



(7) 



yields for the five sites inside the elementary unit centered 
at n: 

x n = a(x n -\ +x n+ i) 

+C {x n -2 + X n+2 ) + P2 X n ±i, 



x n ±i = b (x n + x n ± 2 ) +pi x nT i, 

x n ±2 = a (x„±i + x n ± 3 ) 

+c (x n + x n±4 ) + p 2 x m± , 



(8) 



where we have absorbed the parameters p and z into gen- 
eral "hoping rates" that are initially = = |(1— p), 
c(°) = 0, and p^ = p^ = zp. 

The RG update step consist of eliminating from these 
five equations those two that refer to an odd index, nil. 
After some algebra, we obtain 

X n = b' (x„_ 2 + X n+2 ) + Pi X n ±4, 

X n ±2 = a' (x n + X n±i ) (9) 
+c' (x nT 2 + x n ± 6 ) + p' 2 x m± , 



with 



a = 



b' = 



[aft + c(l-pi)](l+pi) 
\-p\- 2ab 

ab + c(l — p\) 



Pi 



P2 



I -Pi 


- 2ab 


abp 


L 


i-pI- 


-2ab' 


P2(l- 


Pi) 


1 -Pi - 


- 2ab' 


P2(l" 


Pi) 



(10) 



I -p{- 2ab 



If for all sites I = n, n ± 2, n ± 4, . . . in Eq. (9) we further 



identify 1 xi 



C x'^ 2 , we note that the primed equations 



coincide with the unprimed ones in Eqs. (8). Hence, the 
RG recursion equations in (10) are exact at any step k 
of the RG, where unprimed quantities refer to the fcth 
recursion and primed ones to k + 1 . 

Solving Eqs. (10) algebraically at infinite time [which 
corresponds to the limit z / 1, see Eq. (7)] and for 
k+ 1 <~ k — > oo (by dropping the prime on all left-hand pa- 
rameters), we - apparently - obtain only two fixed points 
at a = b = 1/2 and c = p\ = p 2 = 0, and a = b = c = 
and pi = P2 = 1. The first fixed point corresponds to 
an ordinary Id walk without long-range jumps, in the sec- 
ond there is no hopping along the ld-backbone at all and 
the walker stays confined, jumping back-and-forth within 
a single, long-range jump. Yet, both fixed points are un- 
stable with respect to small perturbations in the initial 
parameters. 

Starting with any positive probability p for long-range 
jumps, those dominate over the Id walk at long times. 
Paradoxically, exclusive long-range jumps found at the 
2nd fixed point lead to confinement, itself undermined by 
any positive probability to escape along the Id— line, al- 
lowing to reach even longer jumps. Instead, the process 
gets attracted to a third, stable fixed point hidden inside a 



1 As we will show elsewhere [20], the constant C is determined 
when initial and boundary conditions arc considered, as is essential 
for the case of first transit and return times [21] . 
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Fig. 5: Plot of the probability Pf(A£) of first returns to the 
origin after At update steps on a system of unlimited size. Data 
was collected for three different walks on HN3 with p = 0.1 
(circles), p — 0.3 (squares), and p — 0.8 (diamonds). The data 
with the smallest and largest p exhibit strong transient effects. 
The exact result in Eq. (14), /j, — 1.234 . . ., is indicated by the 
dashed line. 



Fig. 3: Depiction of the (exact) RG step for random walks on 
HN3. Hopping rates from one site to another along a link are 
labeled at the originating site. The RG step consists of tracing 
out odd-labeled variables x n ±i in the top graph and expressing 
the renormalized rates (a', b' , c' , p'i, P2) on the right in terms of 
the previous ones (a, b, c,p\,p2) on the bottom. The node x n , 
bridged by a (dotted) link between x n -i and x n +i, is special 
as it must have n = 2(2j + 1) and is to be decimated at the 
following RG step, justifying the designation of p\. Note that 
the original graph in Fig. 1 does not have the green, dashed 
links with hopping rates (c, c'), which emerge during the RG 



recursion. 
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Fig. 4: Plot of the results from simulations of the mean-square 
displacement of random walks on HN3 displayed in Fig. 1. 
More than 10 7 walks were evolved up to t max = 10 6 steps to 
measure (r 2 )t- The data is extrapolated according to Eq. (2), 
such that the intercept on the vertical axis determines d w 
asymptotically. The exact result from Eq. (13) is indicated 
by the arrow. 



singular boundary layer [17] in the renormalization group 
equations (10) near the confined state. 

We have to account for the asymptotic boundary layer 
in Eqs. (10) with the Ansatz y ~ A y oT k — > for 
y G {a, 6, c, 1 — p\. 1 —P2}, where k — > 00 refers to the fcth 
RG step. Choosing A a = 1, the other A y 's and the eigen- 
values a are determined self-consistently. The only eigen- 
value satisfying the requirement a > 1 is a = 2/cj). Here, 
(j> = (y/5 + l) /2 = 1.6180 ... is the legendary "golden ra- 
tio" [24] defined by Euclid [25]. Hence, every renormaliza- 
tion of network size, L — > L' = 2L, has to be matched by 
a rescaling of hopping rates with a = 2/<f> to keep motion 
along the ld-backbone finite and prevent confinement. 

Extending the analysis to include finite-time corrections 
(i. e., 1 — z <C 1), we extend the above Ansatz to 



y 



(fe) 



A,, 



{\ + {l-z)B y p k + ...} (11) 

for all y S {a, b, c, 1—px, 1— ^2}- In addition to the leading- 
order constants A y and a, also the next-leading constants 
are determined self-consistently, and we extract uniquely 
f3 = 2a. Accordingly, time re-scales now as 



T 



T' = 2aT, 



and we obtain from Eq. (2) with T ~ 
exponent for HN3 

al w = 2 -log 2 = 1.30576. 



(12) 

for the diffusion 
(13) 



The result for d w is in excellent agreement with our sim- 
ulations, as shown in Fig. 4. 

Using the methods from Ref. [21], a far more extensive 
treatment shows [20] that the exponent fi for the proba- 
bility distribution, Pp(At) ~ At _AI , of first-return times 
At is given by 

(j, = 2- — = 1.2342.... (14) 
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The relation between fi and d w is typical also for Levy 
flights [3] , and the result is again borne out by our simula- 
tions, see Fig. 5. It is remarkable, though, that the more 
detailed analysis in Ref. [20] also shows that walks on HN3 
are not uniformly recurrent, as the result of d w > df = 1 
here would indicate. That calculation shows that only 
sites on the highest level of the hierarchy are recurrent. 
While all other sites do share the same exponent \i in 
Eq. (14) for actual recurrences, they have a diminishing 
return probability with decreasing levels in the infinite sys- 
tem limit. This is clearly a consequence of walkers being 
nearly-confined to the highest levels of the hierarchy at 
long times, as expressed by the boundary layer. 

We finally contrast the behavior of HN4 discovered 
above with the analysis of HN3. Clearly, when long- 
range jumps are interconnected as in HN4, there is no 
confinement, the boundary layer disappears [which would 
be similar to a = 1 in Eqs. (11-12) for HN3], and diffu- 
sion spreads ballistically, d w = 1. Our numerical studies, 
and the similarity to Weierstrass random walks [22], fur- 
ther supports that fj, for walks on HN4 is also given by 
Eq. (14), leading to fi = 1. This scaling is again indicative 
of a marginally recurrent state and requires logarithmic 
corrections for proper normalization, as was observed in 
simulations [15]. 



Conclusions. — We conclude with two further con- 
siderations. First, in reference to the potential of these 
networks to interpolate between long-range and a finite- 
dimensional behavior that we invoked in the introduc- 
tion, we just add the following illustrative remark: If the 
probability to undertake a long-distance jump would be 
distance-dependent in each level of the hierarchy, we can 
obtain immediately a new result for walks on HN4 in the 
annealed approximation above. Let p vary with a power 
of the backbone-distance between sites, say p oc r~ a , then 
for each level i of the hierarchy it is r = ri = 2 l , i. e. 
p = pi oc 2~ 4CT , and the weight to make a jump of length 
2 1 in Eq. (5) is given by a = 2 1+CT , leading to d w = 1 + a. 
As can be expected, the analysis of the Weierstrass walk 
breaks down for a — ► 1 _ , at which point the long-range 
jumps become irrelevant and we obtain the results for or- 
dinary Id diffusion. Hence, < a < 1 interpolates analyt- 
ically between long-range and one-dimensional behavior of 
the random walk on HN4. (In fact, the analysis formally 
can be extended to > a > — 1, where the walk becomes 
non-recurrent and is dominated by high levels in the hi- 
erarchy. Yet, the annealed approximation that assumes 
free transitions between different levels of the hierarchy is 
bound to fail.) 

Our final consideration concerns a recent proposal by 
Condamin et al [5] for a very general scaling form for 
mean first-passage times (T) for walks as a function of 
distance r between starting and target site on a graph 
(lattice, network, etc.) of N sites. Based on d w — df, 
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Fig. 6: Plot of the mean first-passage times (T) for walks as a 
function of distance r between starting and target site for HN3 
(top) and HN4 (bottom). The data has been scaled according 
to Eq. (15 such that the data collapses asymptotically onto a 
line that only depends on the model but that is independent 
of system size N. This collapse is excellent for HN3, it is 
somewhat weaker for HN4. Although all system sizes lead to 
linear forms in ln(r), their slope apparently varies with N. This 
could be caused by logarithmic scaling corrections to the slope, 
or by the lack of asymptotic behavior at the available system 
sizes N. 



Ref. [5] determined that 



A + Br d -- d f, d w >d f , 
(T) ~ N {A + Bln(r), d w =d f , 

d w < df, 



(15) 



A - Br d ™- d f 



for fixed constants A, B, independent of N and r. Our 
networks provide a non-trivial set of exponents to explore 
these relations with simple simulations. In particular, 
HN3 with d w — df = 0.30576 provides an instance for 
a powerlaw-divergent mean first-passage time, while HN4 
exactly probes the marginal case d w = dt(= 1) with a 
logarithmic divergence of (T). When plotting (T)/N as 



a function of 



or ln(r), resp., in Figs. 6 we indeed 



obtain a universal straight line over many orders of mag- 
nitude in N and r, indicative of fixed A, B. 
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